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ABSTRACT 


A  static  pressure  probe  was  tested  to  determine  the  feasibility  of 
using  the  probe,  as  an  integral  part  of  a  missile  nose,  to  sense  missile 
altitude.  Experiments  were  conducted  at  Mach  2.0  and  at  Mach  1.51. 

At  Mach  2.0,  the  static  pressure  probe  will  perform  within  altitude 
specifications  of  25,000  feet  ±  2,000  feet  at  angles  of  attack  ranging 
from  -  8  to  ♦  8  degrees.  At  Mach  2.0,  within  an  angle-of-attack  ranging 
from  0  to  6  degrees,  the  probe  will  measure  free  stream  static  pressure 
within  4  percent;  a  4  percent  error  in  measurement  is  equivalent  to  an 
altitude  error  of  900  feet.  The  missile  nose  shock  will  remain  downstream 
of  the  probe  pressure  ports  for  flight  Mach  numbers  above  1.5. 
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I.  INTRODUCTION 


A.  BACKGROUND 

The  ZEPPO  Rocket  is  an  air-launched  system  designed  to  deliver  an 
expendable  pulse  power  communications  jammer  at  long  range.  To  achieve 
this  objective,  the  system  is  loft-launched  at  angles  up  to  45  degrees 
and  at  high  altitude.  At  the  proper  point  in  the  trajectory,  a  deploy¬ 
ment  sequence  is  initiated  by  the  fuze.  At  the  end  of  this  sequence,  the 
expendable  jammer  is  suspended  and  operated  from  a  17-foot  diameter  para¬ 
chute  with  a  descent  rate  of  10  to  12  feet  per  second.  The  ZEPPO  system 
is  5  inches  in  diameter  and  115  inches  in  length,  and  weighs  approximately 
135  pounds.  The  deployment  system  consists  of:  (1)  a  timer  fuze  in  the 
nose;  ( 2 )  an  expulsion  charge;  (3)  a  drogue  parachute ;  (4)  a  main  para¬ 
chute;  and  (5)  an  expendable  electronic  jammer. 

The  rocket  motor  is  the  production  MK71  MOD  1  Zuni  motor.  The  motor 
weighs  80  pounds,  and  is  73.3  inches  in  length  and  5  inches  in  diameter. 

It  has  a  thrust  level  of  6500  pounds  force  for  1.5  seconds,  with  a  burn¬ 
out  spin  rate  of  30  to  33  revolutions  per  second. 

The  fuze,  at  the  present  time,  is  a  modified  MK375  MOD  0  timer  fuze 
with  a  timing  capability  up  to  140  seconds.  It  contains  a  black  powder  ex¬ 
pulsion  charge. 

The  drogue  parachute  has  a  cruciform  geometry  and  is  packed  in  a 
small  bag  that  opens  after  motor  separation.  The  main  parachute  is  a 
lightweight  17-foot  parachute  adapted  from  the  LUU-2/B  flare  system  which 
is  a  standard  production  item.  The  chute  is  packed  in  a  deployment  sleeve 


measuring  5  inches  in  diameter  by  eight  feet  in  length;  this  sleeve  dras¬ 
tically  reduces  snatch  and  opening  loads ,  peimitting  the  lightweight  para 
chute  to  survive  high  "Q"  deployment  without  damage. 

The  optimum  point  for  deployment  initiation  is  25,000  feet,  plus  or 

minus  2,000  feet.  Deployment  at  this  altitude  allows  maximum  payload  time 

aloft,  and  therefore  maximum  operation  time.  System  testing  by  the  Naval 

Weapons  Center,  China  Lake,  California,  has  demonstrated  that  there  is 
0 

difficulty  in  attaining  the  optimum  deployment  altitude  with  the  present 
fuzing  system. 

To  remedy  this  situation,  several  alternatives  were  considered,  one 
of  which  was  a  fuzing  system  initiated  by  a  pressure-sensing  switch,  or 
Baro-Fuze.  Two  alternatives  were  considered  for  measuring  pressure:  a 
shock-swallowing  air  sensor  and  a  static-pressure  probe.  The  Naval  Post¬ 
graduate  School  was  selected  to  evaluate  the  effectiveness  of  the  static- 
pressure  probe.  The  probe  selected  for  this  evaluation  was  developed  by 
Pinckney  [1] . 

Figure  1-1  shows  the  probe  interior  to  the  missile  nose,  where  it 
remains  until  after  launch.  After  rocket  motor  burnout,  and  with  the  mis¬ 
sile  flying  at  small  angles-of-attack,  the  probe  is  extended.  With  the 
probe  extended,  as  shown  in  Figure  1-2,  the  aerodynamic  drag  is  reduced; 
thus  the  decay  of  missile  velocity  is  decreased.  The  range  of  possible 
Mach  numbers  was  considered  to  be  Mach  1.4  to  Mach  2.0  at  small  angles- 


of-attack. 


B.  THESIS  PURPOSE 

The  purposes  of  this  thesis  are:  (1)  to  evaluate  theoretically  the 
static  pressure  probe  general  designs  developed  by  Pinckney  [1]  for  use 
in  the  Mach  1.4  to  Mach  2.0  range;  (2)  to  develop  pressure  distributions 
along  the  probe,  and  select  the  probe  design  best  suited  to  this  Mach 
number  range;  (3)  to  evaluate  experimentally  the  probe  pressure-sensing 
performance  as  an  integral  part  of  the  missile  nose  (with  probe  extended); 
and  (4)  to  report  on  experimentation  in  the  Naval  Postgraduate  School 
supersonic  wind  tunnel  utilizing  the  existing  Mach  1.4  and  Mach  2.0  noz¬ 
zle  blocks. 


II.  THEORY 


A.  STATIC  PRESSURE  PROBE  DESIGN 

Pinckney  [1]  discusses  the  design,  for  supersonic  flow,  of  a  short 
static  pressure  probe  that  is  relatively  insensitive  to  angle -of -attack. 
The  general  probe  design,  which  is  illustrated  in  Figure  2-1,  consists  of 
a  conical  nose  followed  by  a  parabolic  tangent  ogive  fairing  the  nose  into 
a  second,  truncated  cone.  The  second  cone  nates  with  the  cylindrical  after¬ 
body  of  the  probe.  The  parabolic  formula  used  for  the  tangent  ogive  sec¬ 
tion  was  of  the  form: 


Y 

IT 


IA(*/ 


BtV 


♦  C]' 


(2-1) 


The  value  of  Y/Rn  at  the  junction  between  the  tangent  ogive  and  the 
second  truncated  cone  was  specified.  The  conical  nose  half-angle  (B)  and 
the  second  truncated  cone  half-angle  (ca)  were  also  specified.  Based  on 
this  information,  the  constants  A,  B,  and  C  were  calculated.  The  probe 
nose  half-angle  is  dependent  on  the  free-stream  Mach  number,  M»;  the  cone 
angle  must  be  small  enough  so  that  shock  detachment  does  not  occur.  Static 
pressure  holes  are  located  at  points  where  static  pressure  equals  free- 
stream  static  pressure.  These  points  of  pressure  equality  vary  with  the 
length  and  angle  of  the  conical  nose.  Pinckney  [1]  calculated  the  theor¬ 
etical  pressure  distribution  for  this  general  design  utilizing  a  blunt - 
body  program  coupled  with  a  method  of  characteristics  program.  He  made 
calculations  for  B  *  10  degrees  with  u  varying  between  2  and  4  degrees 


for  M.  ■  2.S  and  4.0  and  zero  angle-of -attack.  Wind  tunnel  tests  were  con¬ 
ducted  utilizing  four  different  probe  types.  The  first  two  probe  types  had 
0  ■  10  degrees  and  id  ■  2  or  3  degrees,  and  were  tested  at  M»  *  2.5  and  4. 

The  second  two  probe  types,  which  had  0-20  degrees  and  id  *  2.0  or  3.5 
degrees,  were  also  tested  at  Ms  ■  2.5  and  4.  Tests  were  conducted  at  attack 
angles  ranging  from  0  to  12  degrees. 

Figure  2-2,  which  is  reproduced  from  Kuethe  and  Chow  [2],  is  a  graph 
of  shock  angle  versus  Mach  number  for  various  cone  angles.  Since  the  thrust 
of  this  research  concerns  Mach  numbers  between  1.4  and  2.0,  Figure  2-2 
shows  that  6  cannot  be  greater  than  approximately  15  degrees,  or  shock  wave 
detachment  will  result  at  M»  *  1.4.  Therefore,  a  conical  nose  half-angle 
of  0  «  10  degrees  was  chosen.  Pinckney's  theoretical  results  for  0  ■  10 
degrees,  which  are  shown  in  Figure  2-3,  indicate  that  the  pressure  distri¬ 
bution  which  most  closely  approximates  that  of  the  free-stream  static  pres¬ 
sure  will  be  obtained  with  id  «  3  degrees  at  M*  •  2.5.  Pinckney's  experimental 
results,  which  are  shown  in  Figure  2-4,  indicate  close  agreement  with  the 
theoretical  calculations.  Figure  2-5  illustrates  Pinckney's  probe  design 
for  0  ■  10  degrees  and  id  ■  3  degrees.  Since  Pinckney's  theoretical  calcu¬ 
lations  were  for  much  higher  Mach  numbers,  it  was  considered  prudent,  for 
research  completeness,  to  obtain  theoretical  values  for  the  static  pressure 
distribution  at  lower  Mach  numbers.  For  consistency,  attempts  were  made  to 
obtain  the  computer  program  used  by  Pinckney  [1] ;  however,  the  program  was 
not  immediately  available.  Therefore,  an  alternative  approach  was  sought. 

The  coefficients  of  the  parabolic  tangent  ogive  sections  for  0  ■  10  degrees 
and  id  ■  3  degrees  were  readily  available  and  were  provided  by  Mr.  Pinckney. 

The  values  of  the  coefficients  are  as  follows:  A  -  -0.009267312,  B  «  8.3676155, 


and  C  ■  -402.88107.  The  coefficients  were  verified,  and  Eq.  2-1  provides 
the  ogive  section  body  shape  used-  in  the  theoretical  work  described  in  the 
next  section. 


B.  STATIC  PRESSURE  DISTRIBUTION 

Gawain  and  Schonberger  [3]  discuss  small  pertuxbation  supersonic  flow 
theory  applied  to  a  body  of  revolution.  The  theory  presented  in  Ref.  3  is 
essentially  that  of  classical  slender-body  small -perturbation  theory,  with 
corrected  boundary  conditions  in  the  sense  that  the  boundary  conditions 
are  applied  at  the  body  surface.  In  more  elementary  theory,  the  boundary 
conditions  are  applied  along  the  axis  of  the  body.  The  use  of  the  exact 
boundary  conditions  yields  more  accurate  solutions  than  those  obtained 
from  more  elementary  theory.  The  theory  of  Ref.  3  was  considered  adequate 
for  the  purposes  of  this  thesis,  as  only  small  angles-of-attack  were  con¬ 
sidered.  A  BASIC  language  computer  program,  suitable  for  use  on  a  micro¬ 
processor,  was  developed  based  on  the  theory  and  calculation  procedures 
in  Ref.  3.  The  program,  and  a  sample  data  run  for  *  2.0  at  angle-of- 
attack  o  ■  0  and  10  degrees,  are  given  in  Appendix  A. 

C.  RESULTS  OF  THEORETICAL  CALCULATIONS 
1.  Theoretical  Calculations 

The  pressure  distribution  was  calculated  for  Mo,  »  1.4,  1.6,  1.8, 
2.0,  and  2.S  at  zero  angle-of-attack  (a).  The  M.  »  2.5  calculations  served 
as  a  comparison  between  Pinckney's  theoretical  and  experimental  results  [1] 
and  the  theory  used  in  this  research.  Figure  2-6  is  a  graph  of  P/P„  versus 
X/Rn  for  M*  ■  2.S,  8  ■  10  degrees,  and  u  »  3  degrees.  Comparing  Figures 
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2-3  and  2-6,  it  can  be  seen  that  the  graph  shapes  are  the  sane;  however, 
the  theoretical  results  depicted  in  Figure  2-6  are  slightly  higher  than 
those  of  Figure  2-3  (approximately  2.7  percent).  Yet,  in  comparing  Figure 
2-6  with  the  experimental  results  of  Figure  2-4,  only  a  1  percent  differ¬ 
ence  exists. 

In  addition  to  the  calculations  above,  computations  were  made  for 
angles-of-attack  a  *  1,  2,  4,  6,  8,  and  10  degrees  at  M»  »  1.4,  1.6,  1.8, 
and  2.0.  The  calculation  procedure  and  results  will  be  discussed  in  Sec¬ 
tion  4,  below. 

2.  Probe  Pressure  Hole  Location 

As  was  discussed  previously,  in  Part  A  of  this  chapter,  optimum 
probe  pressure  hole  location  depends  on  the  free-stream  Mach  number.  Due 
to  the  range  of  Mach  numbers,  a  design  compromise  is  made  that  minimizes 
the  pressure  error  over  the  range  of  M*  rather  than  for  one  specific  M». 
Figure  2-7  is  a  graph  of  P/P„  versus  X/R,,  for  ■  1.4,  1.6,  1.8,  and  2.0, 
at  zero  angle-of-attack .  The  graph  shows  clearly  that  the  minimum  differ¬ 
ence  in  P/P.  for  the  range  of  M»  occurs  at  approximately  X/I^,  ■  260.  The 
difference  in  P/P„  due  to  variable  is  approximately  4  percent. 

3.  Error  in  Altitude  Due  to  Error  in  Pressure  Measurement 

Pressure  as  a  function  of  altitude  is  given  by: 

P  «  ?0*<h  .  (2-2) 

where  £  ■  (7.6  km)*1  and  PQ  is  the  pressure  at  sea  level.  Eq.  2-2  was  ob¬ 
tained  by  graphing  pressure  versus  altitude  using  data  from  the  U.  S. 
Standard  Atmosphere,  1978  [4],  and  fitting  a  curve  to  the  data.  Taking  the 
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natural  logarithm  of  both  sides  of  Eq.  2-2  yields: 

In P  ■  -Eh  ♦  ln?Q  .  (2-3) 

Since  PQ  ■  1  atmosphere,  Eq.  2-3  becomes: 

ZnP  -  -Eh  .  (2-4) 

To  find  the  error  in  altitude  due  to  the  difference  in  pressure,  the  deri 
vative  of  Eq.  2-4  is  taken,  which  gives: 

dP/P  -  -Edh  .  (2-5) 


Rearranging  Eq.  2-5,  the  altitude  error  becomes: 


dh  -  -  •  dP/P 


(2-6) 


With  1/E  *  7.62  km  and  dP/P  »  0.04,  the  value  of  dh  is  300  meters,  or  984 
feet.  This  error  is  well  within  the  tolerance  specified  for  probe  perform¬ 
ance,  i.e.,  25,000  feet  ±  2,000  feet. 

4.  Probe  Pressure  Readings  at  Angle-of-Attack 

Figure  2-8  shows  a  cross  section  of  the  static  pressure  probe. 

The  static  pressure  hole  on  the  windward  side  is  number  1,  while  the  hole 
on  the  leeward  side  is  number  3.  Using  the  principle  of  conservation  of 
mass,  the  mass  flow  rate  in,  m,,  must  equal  the  mass  flow  rate  ouc,  m  ; 

S'  O 

therefore: 


f 
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PjVjAj  +  2p2V2A2  ♦  p3V3A3  »  0  .  (2-7) 

Since  Pj  *  P2  *  P3  and  Ai  *  A2  *  A3*  E<1*  2"7  r®duces  t0: 

V3  +  2V2  +  V3  -  0  .  (2-8) 

Note  that  Eqs.  2-7  and  2-8  require  the  velocity  to  be  an  algebraic  quan¬ 
tity:  a  flow  into  the  probe  is  defined  as  positive.  Streeter  and  Wylie  [S] 
demonstrate  that  flow  through  an  orifice  can  be  represented  by: 

P0  -  P  *  *SCPV2  ,  (2-9) 

where  PQ  is  the  pressure  upstream  of  the  orifice,  P  is  the  pressure  down¬ 
stream  of  the  orifice,  and  C  is  a  constant  for  incompressible  flow.  Writing 
Eq.  2-9  in  terms  of  the  parameters  of  Figure  2-8  yields: 


I»1  -  Ppl  • 

(2-10) 

|p,  -  Ppl  •  >»2v£2 

(2-11) 

|P3  -  Ppl  •  >»3v*c,  , 

(2-12) 

where  Pp  is  the  pressure  measured  by  the  probe  and  C  is  the  flow  coeffici¬ 
ent  for  pressure  drop  across  the  pressure  holes.  Solving  Eqs.  2-10,  2-11, 
and  2-12  for  Vj,  V2,  and  V3,  respectively,  and  substituting  into  Eq.  2-8, 
yields: 

[2\PrPp\]h  i  2[2|P2-Pp|]‘s  ±  [2\P3-Pp\]H  -  0  .  (2-13) 

Simplification  to  this  form  was  made  possible  because  P^  s  P2  *  P3  and 
ci  m  c2  a  c3*  Note  that  the  value  of  C  is  not  needed,  since  C  cancels  in 
Eq.  2-13. 
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Eq.  2-13  is  nonlinear,  and  its  eight  mathematically  possible  cases 
are  shown  in  Table  2.1.  A  plus  sign  (+)  indicates  flow  into  the  probe;  a 
minus  sign  (-)  indicates  flow  out  of  the  probe.  Some  cases  which  are  mathe¬ 
matically  possible  but  not  physically  possible  can  be  eliminated. 


Case 

P.-P 

P  -P„ 

"O 

i 

•o 

1  P 

2  P 

3  p 

1 

+ 

+ 

+ 

2 

- 

+ 

+ 

3 

+ 

- 

+ 

4 

+ 

- 

5 

+ 

- 

- 

6 

- 

+ 

- 

7 

- 

- 

+ 

8 

_ 

- 

Table  2.1  Mathematically  Possible  Cases  for  Eq.  2-13 

Cases  1  and  8  can  be  eliminated  immediately  as  not  satisfying  con¬ 
servation  of  mass.  In  addition,  at  positive  angles-of-attack  (i.e.,  probe 
tip  up),  the  pressure  Pj  will  be  higher  than  ?2  or  P3»  Therefore,  cases  2, 
6,  and  7  may  be  ruled  out  as  not  possible.  An  iterative  computer  program 
was  written  to  solve  Eq.  2-13  for  P^,  given  that  cases  3,  4,  and  5  are 
the  only  possible  cases.  A  listing  of  the  program  and  a  sample  data  run 
can  be  found  in  Appendix  B. 

Figure  2-9,  which  plots  P/Pa  as  a  function  of  a,  shows  probe  theo¬ 
retical  performance  for  values  of  a  ranging  between  zero  and  ten  degrees. 
The  data  in  Figure  2-9  were  calculated  for  *  1.4  and  for  M*  ■  2.0. 
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III.  EXPERIMENTAL  APPARATUS 


WIND  TUNNEL 


1.  Mach  2.0  Wind  Tunnel  Tests 


Wind  tunnel  tests  were  conducted  in  the  Naval  Postgraduate  School 
4x4  inch  supersonic  wind  tunnel  using  the  existing  Mach  2.0  nozzle  blocks 
(The  Mach  number  of  2  is  nominal;  the  actual  Mach  number  is  1.92).  Figures 
3-1  and  3-2  are  two  photographs  of  the  general  experimental  apparatus.  They 
present  two  views  of  the  supersonic  wind  tunnel  and  its  associated  instru¬ 
mentation.  Utilizing  a  one-fifth  scale  model,  and  regulating  wind  tunnel 
stagnation  pressure,  an  actual  Reynolds  number  of  Re  *  1.3  x  10s  was  at¬ 
tainable,  which  compares  favorably  with  the  Reynolds  number  for  the  full- 
scale  probe  at  25,000  feet  (viz..  Re  ■  4.3  X  10s).  The  reference  length 
for  the  calculation  of  the  Reynolds  number  is  0.020  meters.  Twenty-three 
wind  tunnel  tests  were  conducted  at  M»  *  1*92  with  a  ■  0,  1,  2,  4,  6,  8, 
and  10  degrees. 

Prior  to  any  testing,  the  fine-mesh  wire  screen  in  the  wind  tun¬ 
nel  plenum  was  discovered  to  be  tom,  and  hanging  loose  in  the  plenum 
chamber.  The  screen  serves  to  insure  a  flat  velocity  profile  (flat  stag¬ 
nation  pressure  profile)  and  to  generate  small-scale  turbulence  which  dis¬ 
sipates  prior  to  arrival  at  the  test  section.  The  plenum  chamber  exterior 
can  be  seen  in  Figure  3-1.  Sufficient  manpower  to  replace  the  plenum  cham¬ 
ber  screen  was  not  available.  Therefore,  the  screen  was  removed .  To  verify 
that  the  flow  in  the  wind  tunnel  test  section  was  uniform,  an  additional 
pressure  tap  was  installed  opposite  one  of  the  existing  pressure  taps. 


Testing  showed  no  significant  pressure  difference  (approximately  0.5  per¬ 
cent)  between  the  two  pressure  taps.  Therefore,  the  pressure  measurements 
suggest  that  the  test  section  velocity  profile  was  reasonably  uniform  at 
Mach  2.0. 

2.  Mach  1.4  Wind  Tunnel  Tests 

Tests  were  attempted  utilizing  the  existing  Mach  1.4  nozzle  blocks. 
Supersonic  flow  could  not  be  established  in  the  test  section.  Appendix  C 
provides  a  sequential  listing  of  testing  attempts  and  corrective  measures 
taken  between  tests,  and  discusses  the  three  most  probable  reasons  for 
testing  failure. 

3.  Mach  1.51  Wind  Tunnel  Tests 

Since  experimental  data  at  a  lower  Mach  number  were  still  desired, 
various  alternatives  were  explored.  The  four  alternatives  investigated  were: 
(1)  to  replace  the  plenum  chamber  screen;  (2)  to  design  and  manufacture 
Mach  1.5  nozzle  blocks;  (3)  to  enlarge  the  existing  Mach  1.4  nozzle  block 
second  throat  by  machining;  or  (4)  to  attempts  tests  using  a  combination 
of  Mach  1.4  and  Mach  2.0  nozzle  blocks.  Alternatives  1  through  3  were  ruled 
out  due  to  inadequate  manpower,  funding,  and/or  research  time.  Alternative 
4  was  attempted.  Testing  in  this  unorthodox  manner  was  attempted  in  order 
to  verify  that  the  missile  nose  shock  did  not  detach  at  a  lower  Mach  num¬ 
ber.  Two  shock  waves  occur  at  the  nose  of  the  missile  with  the  probe  extended. 
A  nearly  conical  shock  is  attached  to  the  probe.  A  curved  shock  wave  origi¬ 
nates  due  to  the  missile  ogive;  this  curved  shock  wave  is  termed  the  missile 
nose  shock.  The  probe  conical  shock  and  missile  nose  shock  merge  to  form 
the  bow  shock  wave.  The  probe  conical  shock,  missile  nose  shock,  and  re¬ 
sulting  bow  shock  wave  can  be  seen  in  Figure  3-3,  which  is  a  photograph  of 


the  model  mounted  in  the  wind  tunnel.  The  probe  is  oriented  at  zero  angle- 


of-attack,  and  the  flow  is  at  »  1.92.  Both  the  probe  conical  shock  wave 
and  the  missile  nose  shock  wave  are  attached.  If  the  missile  nose  shock 
moves  upstream  of  the  probe  pressure  ports,  erroneous  pressure  readings 
are  obtained. 

The  Mach  number  gradient  in  the  wind  tunnel  test  section  resulting 
from  the  combination  of  nozzle  blocks  was  calculated  using  the  method  of 
characteristics  found  in  Liepmann  and  Roshko  [6] .  Figure  3-4  shows  the 
Mach  number  gradient  in  the  test  section  and  the  shock  waves  that  result 
using  the  combination  of  nozzle  blocks.  Figure  3-5  is  a  photograph  of  the 
Mach  2.0  and  Mach  1.4  nozzle  blocks  installed  in  the  wind  tunnel.  The  theo 
retical  Mach  number  obtained  using  this  combination  of  nozzle  blocks  was 
M,,,  =  1.51  at  the  probe  nose,  *  1.538  at  the  top,  and  M^  *  1.16  at  the 
bottom  of  the  wind  tunnel  test  section.  Experimental  results,  which  are 
discussed  in  Chapter  IV,  were  used  to  verify  theory. 

B.  MODEL  DESIGN 

A  one-fifth  scale  model  of  the  static  pressure  probe  and  missile  nose 
was  used.  Since  a  complete  model  could  not  be  mounted  in  the  tunnel,  the 
missile  nose  was  truncated  at  40  percent  of  the  nose  length.  Maximum  model 
size  was  dictated  by  maximum  allowed  wind  tunnel  blockage  (approximately 
6.4  percent)  at  Mach  1.4.  The  allowed  wind  tunnel  blockage  limited  the 
combined  frontal  area  of  the  model  and  mounting  system.  Pinckney  [1]  sup¬ 
plied  the  drawings  for  the  static  pressure  probe.  The  Naval  Weapons  Center 
China  Lake,  California,  provided  the  formula  and  coefficients  for  the  Von 
Karmen  ogive  missile  nose.  The  Von  Karmen  ogive  is  given  by: 
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(3-1) 


Y  *  ■£=  •  [cos”1  (1  -  -  Jjsinwcos^U  -  Q))]h 

for  the  full-scale  missile,  where  R  *  2.5  inches  and  L  =  10  inches.  Figure 
3-6  is  a  front  view  of  the  model  and  double-wedge  airfoil  used  to  mount 
the  model  in  the  wind  tunnel.  Design  compromises  between  the  maximum  mis¬ 
sile  nose  diameter  and  maximum  wedge  thickness  were  required  to  remain 
within  wind  tunnel  blockage  limits.  The  maximum  diameter  of  the  model  nose 
used  is  0.792  inches,  and  the  maximum  wedge  thickness  is  0.165  inches.  The 
mounting  wedge  length  is  3.900  inches,  and  a  wedge  half-angle  of  9  degrees 
was  utilized  to  prevent  oblique  shock  detachment  at  Mach  1.4.  The  frontal 
area  of  the  model  plus  mounting  system  was  1.005  square  inches,  which  is 
within  the  allowed  wind  tunnel  blockage.  Figure  3-7  is  a  photograph  of 
the  model  ogive  nose  and  the  extended  static  pressure  probe.  The  ogive  is 
1.200  inches  in  length.  The  pressure  probe  extends  0.680  inches  from  the 
missile  nose.  Figure  3-8  is  a  photograph  of  the  static  pressure  probe 
mated  to  the  missile  nose.  The  four  probe  pressure  holes  are  90  degrees 
apart  and  0.0135  inches  in  diameter,  and  are  located  0.1325  inches  from 
the  probe  conical  nose  tip.  The  parabolic  tangent  ogive  and  second  trun¬ 
cated  cone  start  at  0.0525  inch  and  0.1225  inch,  respectively,  from  the 
probe  nose  tip,  and  are  made  of  0.060-inch  outside  diameter,  0.040-inch 
inside  diameter  stainless  steel  tubing.  The  probe  nose  cone  through  the 
second  truncated  cone  is  made  of  0.0625-inch  diameter  hole  on  centerline. 
After  machining,  the  drill  rod  is  mated  to  the  stainless  steel  tubing 
with  silver  solder.  Figure  3-9  is  a  photograph  of  the  top  view  of  the 
model  and  mounting  wedge.  Figure  3-10  is  a  photograph  of  the  model  mounted 
in  the  wind  tunnel,  and  Figure  3-11  is  a  close-up  of  the  model  mounted  in 
the  wind  tunnel  test  section. 
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C.  INSTRUMENTATION 


Pressure  taps  in  the  wind  tunnel  wall  and  on  the  model  surface  were 
connected  to  a  Giannini  12-position  scanivalve  with  Tygon  tubing.  The 
scanivalve  has  a  range  from  absolute  zero  to  100  PSIA.  The  scanivalve  was 
connected  to  a  Gould  Statham  absolute  pressure  transducer.  Model  PA732TC- 
100-350.  The  transducer  was  powered  by  a  five-volt  excitation;  it  has  a 
pressure  range  from  absolute  zero  to  100  PSIA,  with  a  maximum  frequency 
response  of  8,700  Hz.  A  variable-gain  differential  amplifier  conditioned 
the  transducer  output  to  an  Analog  Devices  4*s-digit  voltmeter.  The  voltage 
displayed  is  directly  proportional  to  the  pressure.  Transducer  calibration 
was  accomplished  utilizing  a  60-inch  mercury  column  and  vacuum  pump.  Por 
higher  pressure,  a  Wallace  and  Tiernan  bourdon  tube  absolute  pressure  gauge 
was  used  (range  absolute  zero  to  100  PSIA) .  Schlieren  photographs  of  the 
ow,  including  shock  waves,  were  taken.  A  1,000-watt  mercury  vapor  light 
source  and  a  reflex-type  camera  converted  to  accept  Polaroid  Type  152  film 


were  used. 


IV.  EXPERIMENTAL  RESULTS 


A.  MACH  2.0  EXPERIMENTAL  RESULTS 
1.  Data  Reduction  Technique 

Twenty- three  wind  tunnel  tests  were  conducted  at  Mach  2.0  with 
angles-of-attack  varying  from  -  2  to  +  10  degrees.  Tables  4.1  and  4.2 
depict  sample  wind  tunnel  test  data  for  a  *  0  and  2  degrees,  respectively. 

Run:  4  P  *  29.85  in  Hg  Date:  24  September  1980 
^  *  1.92  a  *  0°  T  *  16.2°C  Time:  0902 


Reading 

Model 

(PSIG) 

Port  n 
(PSIG) 

Pieman 

(PSIG) 

1 

-9.37 

-9.25 

+23.12 

2 

-9.37 

-9.30 

+22.82 

3 

-9.63 

-9.49 

+22.14 

4 

-9.77 

-9.77 

+19.92 

5 

-9.86 

-9.82 

+19.27 

Table  4.1  Sample  Wind  Tunnel  Data  for  a  »  0  Degrees  and  M»  *  1.92 


Run:  6  P  »  29.88  in  Hg  Date:  24  September  1980 
Mo,  *  1.92  a  «  2°  T  ■  20.5°C  Time:  1027 


Reading 

Model 

(PSIG) 

Port  n 
(PSIG) 

Plenum 

(PSIG) 

1 

-9.12 

-9.02 

+24.75 

2 

-9.16 

-9.06 

+24.58 

3 

-9.18 

-9.10 

♦24.37 

4 

-9.21 

-9.18 

+24.14 

S 

-9.35 

-9.28 

♦23.90 

Table  4.2  Sample  Wind  Tunnel  Data  for  a  *  2  Degrees  and  M»  ■  1.92 


The  ratio  of  measured  static  probe  pressure  to  wind  tunnel  wall 


pressure,  P/P.,  was  calculated  for  each  set  of  data  points.  The  arithmetic 
mean,  variance  (N  weighting),  and  standard  deviation  (N  weighting)  were 
calculated  using  the  calculation  procedure  found  in  Haber  and  Runyon  [7] . 
Table  4.3  gives  these  values  for  various  angles-of-attack,  a. 


Angle-o f- 
Attack 

Theoretical 

P/P. 

Experimental 

P/P. 

Experimental 

Standard 

Deviation 

6* 

1.003 

•  0.987 

±0.00861 

1° 

1.001S 

0.984 

±0.0102 

2° 

0.996 

0.987 

±0.0097 

4° 

0.976 

0.978 

±0.0122 

6° 

0.952 

0.974 

±0.008 

8° 

0.938 

0.932 

±0.013 

© 

O 

H 

0.918 

0.927 

±0.016 

Table  4.3  Angles-of -Attack  Versus  Theoretical  and 
Experimental  P/P.  for  M .  »  2. 0 

2.  Static  Pressure  Probe  Performance 

Figure  4-1  is  a  plot  of  the  data  listed  in  Table  4.3.  Calculated 
points  are  indicated  by  the  letter  X,  while  experimental  mean  values  are 
denoted  by  an  asterisk  (*)  and  standard  deviations  are  denoted  by  the 
character  The  scale  of  the  ordinate  emphasizes  greatly  any  deviation 
from  P/P.  *  l.  Note  that  an  experimental  error  bar  is  not  shown  for  angle- 
of-attack;  such  an  error  bar  would  be  horizontal.  Pressure  measurements 
for  angles-of-attack  of  -  1  and  -  2  degrees  are  included  in  the  measure¬ 
ments  shown  for  o  ■  1  and  2  degrees.  A  change  in  a  from  +  2  to  -  2  degrees 
had  little  or  no  effect  on  pressure  measurements.  Figure  4-3  shows  that 


probe  performance  is  reasonably  constant  in  the  range  o  »  0  to  6  degrees; 
the  corresponding  values  of  mean  P/P.  are  0.987  and  0.974.  For  angles -of - 
attack  greater  than  6  degrees,  the  measured  pressure  ratio  is  1.6  percent 
less  than  that  predicted  by  the  theory.  The  mean  value  of  the  measured 
pressure  ratio  is  almost  exactly  that  of  the  theory  at  a  ■  4  degrees;  at 
o  *  6  degrees,  the  measured  P/P.  is  2.6  percent  greater  than  theoretical 
predictions. 

Using  Eq.  2-6  with  dh  >  2,000  feet  (609.6  meters)  and  1/g  >7.62 
kilometers,  the  fractional  change  in  pressure  is  0.08.  Therefore,  taking 
experimental  uncertainty  into  consideration,  the  probe  will  perform  within 
altitude  specifications  for  angles-of-attack  less  than  or  equal  to  8  de¬ 
grees  at  Mach  2.  As  stated  in  Chapter  III,  Figure  3-3  is  a  photograph  of 
the  model  mounted  in  the  wind  tunnel.  The  probe  is  oriented  at  zero  angle- 
of-attack,  and  the  flow  is  at  M.  *  1.92.  Both  the  probe  conical  shock  wave 
and  the  missile  shock  wave  are  visible  and  attached.  The  horizontal  lines 
just  above  the  model  are  due  to  condensation  within  the  boundary  layer  on 
the  inside  of  the  test  section  window.  The  dark  vertical  line  aft  of  and 
below  the  model  is  the  shadow  of  the  Tygon  tubing  used  to  measure  probe 
pressure.  Figure  4-2  is  a  photograph  of  the  model  oriented  at  o  ■  ♦  1 
degree  within  the  wind  tunnel  operating  at  M.  «  1.92. 

B.  MACH  1.51  EXPERIMENTAL  RESULTS 

Seven  wind  tunnel  tests  were  conducted  using  the  combination  of  Mach 
1.4  and  Mach  2.0  nozzle  blocks.  The  model  and  plenum  pressures,  as  well  as 
pressures  at  five  static  pressure  ports  located  along  the  wind  tunnel  wall, 
were  recorded  for  each  wind  tunnel  test.  The  static  pressure  ports  extended 
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from  just  aft  of  the  nozzle  throat  to  the  test  section.  The  position  of 
the  static  pressure  ports  relative  to  the  wind  tunnel  test  section  is 
shown  in  Figure  3-1.  The  static  pressure  measurements  were  used  to  verify 
the  calculations  based  on  the  method  of  characteristics  referred  to  in 
Chapter  III.  The  measured  static  pressure  for  the  seven  wind  tunnel  tests 
was  recorded,  grouped  by  test  point,  and  the  arithmetic  mean  and  standard 
deviation  were  calculated  for  each  point.  Theoretical  values  and  actual 


experimental  values  were 

compared,  and  are 

shown  in  Table  4 

.4. 

Theoretical 

Observed  Mean 

Percent 

Test 

Pressure 

Pressure 

Theoretical 

Pressure 

Point 

(PSIA) 

(PSIA) 

Mach  Number 

Difference 

pi 

P 

2 

5.97 

5.60 

1.57 

5.9 

6.05 

6.35 

1.56 

4.7 

P3 

7.41 

7.92 

1.42 

6.5 

P 

4 

8.87 

9.12 

1.29 

2.7 

P5 

10.00 

10.89 

1.20 

8.1 

Table  4.4  Theoretical  and  Observed  Pressure  for  a 
Combination  of  Mach  1.4  and  Match  2.0  Nozzle  Blocks 


At  Mach  1.5,  a  pressure  difference  of  approximately  7  percent  results 
in  a  Mach  number  change  of  only  0.05.  Therefore,  the  calculated  Mach  num¬ 
bers  are  accurate  within  ±  0.05,  and  a  predicted  Mach  number  of  M*  *  1.50 
±  0.05  at  the  probe  tip  appears  reasonable.  Further  evidence  supporting  a 
value  of  M*  ■  1.5  can  be  seen  in  Figure  4-3,  which  is  a  Schlieren  photo¬ 
graph  of  the  model  in  the  wind  tunnel  test  section  during  one  of  the  test 
runs.  The  Mach  waves  emanating  from  the  top  and  bottom  of  the  test  section 
are  oriented  at  38  and  46  degrees,  respectively.  At  Mach  angle  w  *  46 


degrees,  the  value  of  M.  is  1.38.  The  probe  conical  nose  shock  angles  are 
<^1  ■  38  (top)  and  $2  ■  46  (bottom),  indicating  that  the  probe  is  at  angle- 
of -attack.  A  rough  estimate  of  the  angle-of-attack  is  given  by: 


With  *  38  and  •  46,  the  value  of  a  is  4  degrees.  An  accurate  value 
for  a  would  require  further  study.  An  angle-of-attack  between  4  and  6  de¬ 
grees  was  predicted  theoretically.  Since  all  shocks  in  the  test  section 
are  curved  due  to  Mach  gradient,  all  Mach/shock  wave  angles  were  measured 
locally,  i.e.,  in  close  proximity  to  the  physical  object  generating  them. 
In  addition.  Figure  4-3  shows  a  shock  wave  just  forward  of  the  probe  tip; 
this  shock  wave  was  predicted  theoretically,  and  is  shown  in  Figure  3-4. 

Due  to  the  Mach  number  gradient  and  confused  flow  in  the  test  section 
the  pressure  measured  by  the  probe  is  not  accurate,  and  hence  is  not  rele¬ 
vant.  Table  4.5,  which  presents  sample  wind  tunnel  test  data  with  the 
asymmetric  nozzle  blocks,  is  included  for  completeness  only.  The  unortho¬ 
dox  wind  tunnel  tests  did  prove,  however,  that  the  missile  nose  shock  will 
not  detach  at  Mach  numbers  above  1.5. 

Figure  4-4  is  a  graph  of  distance  versus  theoretical  Mach  number  and 
pressure  (PSIA)  in  the  wind  tunnel  test  section.  The  theoretical  pressure 
is  indicated  by  the  character  M0".  Pressures  measured  by  the  probe,  and 
the  Mach  numbers  calculated  from  the  Mach  waves  emanating  from  the  wind 
tunnel  walls,  are  also  indicated.  Moreover,  Mach  waves  at  several  points 
in  the  plane  of  the  probe  nose  were  measured  and  recorded  on  the  graph. 
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Table  4.5  Sample  Wind  Tunnel  Teat  Data  for  a  =  4  Degrees  and  M„.  -  1.51 


V.  CONCLUSIONS  AND  RECOMMENDATIONS 


A.  CONCLUSIONS 

The  following  conclusions  result  from  the  static  pressure  probe  study: 

(1)  The  static  pressure  probe  of  Ref.  1,  with  0  ■  10  degrees  and 
at  ■  3  degrees,  will,  as  an  integral  part  of  the  aissile  nose, 
perform  within  altitude  specifications  of  25,000  feet  ±  2,000 
feet  at  Mach  2.0  at  angles-of-attack  from  -  8  to  ♦  8  degrees. 

(2)  At  Mach  2.0,  within  an  angle -of- attack  range  froai  0  to  6  degrees, 
the  probe  will  measure  free -stream  static  pressure  within  4  per¬ 
cent.  A  4  percent  error  in  measurement  is  equivalent  to  an  alti¬ 
tude  error  of  900  feet. 

(3)  The  missile  nose  shock  remains  downstream  of  the  probe  pressure 
ports  for  flight  Mach  numbers  above  1.5. 

(4)  Since  the  theory  predicted  probe  performance  within  2.6  percent 
at  Mach  2.0,  and  since  the  theory  predicted  satisfactory  per¬ 
formance  at  Mach  1.4,  the  static  pressure  probe  should  perform 
within  specifications  at  Mach  1.4. 

B.  RECOMMENDATIONS 

The  minimum  operating  Mach  number  for  the  static  pressure  probe  de¬ 
pends  on  the  position  of  the  missile  nose  shock  wave.  At  some  lower  Mach 
number,  the  missile  nose  shock  wave  moves  upstream  of  the  pressure  ports. 

M  is  defined  as  the  minimum  Mach  number  for  satisfactory  probe  operation. 


The  potential  use  of  the  static  pressure  probe  as  an  altitude-measuring 

device  is  limited  by  the  minima  Mach  number,  M  ,  and  by  excursions  in  the 

m 

angle-of-attack.  Recooaended  further  studies  are  as  follows: 

(1)  Replace  the  aissing  supersonic  wind  tunnel  plenua  chamber  screen, 
design  and  manufacture  a  set  of  Mach  1.4  nozzle  blocks  with  a 
larger  second  throat,  and  test  the  probe.  Determine  whether 

M_  <  1.4  or  M  >  1.4. 

7ft  7ft 

(2)  Once  the  actual  flight  Mach  number  has  been  determined,  recalcu¬ 
late  the  theoretical  distribution,  if  required,  and  reposition 
the  probe  pressure  ports  to  remove  the  error  due  to  Mach  number 
range — an  error  presently  inherent  in  probe  design. 


3°  Half  , 
Conical  S 


1-1  Static  Pressure  Probe  Interior  to  Miesile  Nose 


cone  section 


Figure  2-1  General  Static  Pressure  Probe  Design  (Reproduced  from  Pinckney  [lp 


Figure  2-2  Shook  Angle  Versus 
(Reproduced  from 


taoh  Humber  for  Various  Cone  Angles 
l lethe  and  Chow  [2 ]) 
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Theoretical  Static  Pressure  Distribution  as  a  Function  of  X/R^  for  =  2.5, 
1  =  10  Degrees ,  and  u  -  3  Degrees  (Reproduced  from  Pinckney  [1] J 


Figure  2-5  Static  Pressure  Probe  Design  for  3  »  10  Degrees  and  a>  *  3  Degrees  (Reproduced  from  Pinckney  [1]/) 


I  Probe  Static  Pressure  Distribution  for  *  2. 5,  0  *  10  Degrees t  and  u>  =  3  Degrees 


X/RN 

Theoretical  Probe  Static  Pressure  Distribution  for  = 

at  Zero  Angle-of -Attack 


l.ai 


45 


Figure  2-9  Theoretical  Probe  Static  Pressure  Distribution  as  a  Function  of  Angle-of -Attack 
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Instrumentation  (Rear  Quarter  View) 


Instrumentation  (Front  View) 
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Model 


Compression  waves  i 
to  form  an  oblique 


and  2.0  (Top)  Nozzle  Blocks 


■d  the 


•  T  ’■ 


3-10  Photograph  of  the  Model  Mounted  in  the  Wind  Tunnel 
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Figure  3-11  Close-Up  View  of  the  Model  Mounted  in  the  Wind  Tunnel 


Sohlieren  Photograph  of  the  Model  Oriented  at  a  =+  1  Degree  In 
the  Wind  Tunnel  Operating  at  Moo  =  1-92 


Figure  A-l  Sahlieren  Photograph  of  Shook  in  Wind.  Tunnel  Test  Section  During  an  Attempt 
to  Obtain  Supersonic  Flow  with  the  Mach  1.4  Nozzle  Blocks 


APPENDIX  A:  PRESSURE  DISTRIBUTION  PROGRAM  LISTING  AND  SAMPLE  DATA  RUNS 

Note:  In  this  typed  listing,  and  in  the  sample  data  runs,  the  character  " 
signifies  "greater  than,"  and  "±"  signifies  "less  than."  The  ampersand  (&) 
denotes  exponentiation. 

00001  LPRINT"SUPER SONIC  FLOW  PRESSURE  COEFFICIENTS  ABOUT  A  YAWED," 

00002  tPRINT'SLENDER,  JOINTED,  BODY  OF  REVOLUTION.  PROGRAM  BY" 

00003  LPRINT'J.R.SCHONBERGER  AND  K.D.TILLOTSON.  BASED  ON  " 

00004  LPR INT 'CALCULATION  PROCEDURES  BY  PROF.  T.H.GAWAIN." 

00005  LPRINT'NAVAL  POSTGRADUATE  SCHOOL,  MONTEREY  CA.  FALL  1980" 

00006  CLEAR 
00010  CLS:  PI=3.14159 

00020  DIM  AC20,20),BC20,20),CC20,20),CK20),D(20,20),FK20),GC20),GAC2 
0),GB(20),GCC20),GN(20),GOC20),GK20),IX(20),JOC20),JRC20),JXC20 
),RC20),RPC20),TC20,20),XC20),PRC20) 

00025  LPRINT"ENTER  XI"; : IN^TT  XI rLPRINT  XI 

00027  LPRINT"ENTER  X2"; : INPUT  X2:LPRINT  X2 

00029  LPRINT"ENTER  X3"; : INPUT  X3:LPRINT  X3 

00031  tPRINT"ENTER  B1";:IN>UT  B1:LPRINT  B1 

00033  LPRINT'ENTER  B2";rIW»UT  B2 rLPRINT  B2 

00035  LPRINT"ENTER  COEFF  A"; : IM»UT  AA: LPR I NT  AA 

00037  LPRINT'ENTER  COEFF  B"; :  INPUT  BB rLPRINT  BB 

00039  LPRINT’ENTER  COEFF  C";  r  IW>UT  CC:LPRINT  CC 

00041  tPRINT'ENTER  RX2"; : INPUT  RX rLPRINT  RX 

00043  LPRINT'ENTER  FREE  STREAM  MACH  NBR";rIN»UT  M rLPRINT  M 

00045  MU=ATN(  1  /M/SQRC-  l/M54 1  /M+-1 )  ) 

00047  RY=RX+TANCMU)xCX3-X2) 

00049  LPRINT"ENTER  NBR  OF  SUB INTERVALS"; : INPUT  NrLPRINT  N:J=N-1 
00051  LPRINT'ENTER  LENGTH  OF  SUBINTERVAL"; : INPUT  H rLPRINT  H 


00053  FOR  Z=1T0J 

00055  XCDs^H/Cl-TANCBO/TANCMU)) 

00057  R(Z)=X(Z)*TAN(B1) 

00059  RPCZ)=TANCB1) 

00061  IF  XCZ)±=X1  THEN  NEXT  Z  ELSE  63 

00063  FOR  Y=ZTOJ 

00065  R=1-AA/(TAN(MJ))S2 

00067  S=-1S«C2S*Y}«H+BB/CTANCMU))S2) 

00069  T=-1J{CC/(TANCMU))S2+CY}«H)S2 
00071  XCY)=C-S+SQRCSS2-45<RJiT))/C2}tR) 

00073  RCY)=SqRCAA«XCY)S2+BB!tXCY)+CC) 

00075  RPCY)=C2*AA!!X(Y)+BB)/(2KR(Y)) 

00077  IF  XCY)±=X2  THEN  NEXT  Y  ELSE  79 
00079  FOR  Z=YTOJ 

00081  XCZ)=CCRX-X25iTAN<B2))/TANCMU)+Z«H)/Cl-TANCB2)/TANCMU)) 
00083  RCZ)=TANCB2>*SX(Z)+RX-TANCB2)*X2 
00085  RP(Z)=TANCB2) 

00087  IF  X(Z)±=X3  THEN  NEXT  Z  ELSE  89 

00089  FOR  Y=ZTOJ 

00091  X(Y)=RY/TANCMU)+Y*H 

00093  R(Y)=RY 

00095  RPCY)=0 

00097  NEXT  Y 

00099  LPRINT  TABa),M";TABCl0)"X(J)'}TABC20)',RCJ),,;TABC30),,RPCJ)" 
00101  FOR  Z=1T0J 

00103  LPRINT  TABC1)Z;TABC10)XCZ);TABC20)RCZ);TABC30)RPCZ) 
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00110  IF  Z=10  THEN  INPUT”PRESS  ENTER  TO  CONTINUE”; N$ 

00120  NEXT  Z 
00190  B=SQR(MS2-1) 

00192  ON  ERROR  GOTO  1200 
00195  U=0 
00200  FOR  Z=1T0J 
00205  U=U+1 
00210  FOR  Y=0TOU 

00220  T(Z,Y)=SQR((X(Z)-Y}40&2-(BS{R(Z))82) 

00230  NEXT  Y 
00231  NEXT  Z 
00235  U=0 
00240  FOR  Z=1T0J 
00245  U=U+1 
00250  FOR  Y=1T0U 

00260  A(Z/Y)=LOGCCCX(Z)-CY-1)scH)+TCZ/Y-1))/CCXCZ)-Y!JH)+TCZ,Y))) 

00270  NEXT  Y 

00275  NEXT  Z 

00280  FOR  Z=1T0J 

00290  FOR  Y=lTOvJ 

00300  BCZ,Y)=CTCZ/Y-1)-TCZ/Y))/RCZ) 

00310  NEXT  Y 

00311  NEXT  Z 

00320  FOR  Z=1T0J 

00330  FOR  Y=1T0J 

00340  CCZ,Y)=B(Z,Y)-RP(Z>'!A(Z,Y) 


00350  NEXT  Y 
00351  NEXT  Z 
00355  U=0 

00360  FP(1)=RP(D/CC,1) 

00370  FOR  Z=2  TO  J 
00380  U=U+1 
00390  FOR  Y=1T0U 
00400  W=W+C(Z,Y)«FPCY) 

00410  NEXT  Y 

00420  FP(Z)=(RP(Z)-VO/C(Z,Z) 

00430  W=0 
00440  NEXT  Z 
00450  FOR  Z=lTOd 
00460  FOR  Y=1T0J 
00470  W=W+A(Z,Y)ieFP(Y) 

00480  NEXT  Y 
00490  IX(Z)=W 
00500  W=0 
00510  NEXT  Z 
00520  FOR  Z=1T0J 
00540  NEXT 
00550  FOR  Z=1T0J 

00560  GO(Z)=a-Cl+RPCZ)82)«(l-IX(Z))S2+CMXlX(Z))S2) 
00565  NEXT  Z 
00610  FOR  Z=1  TO  J 
00620  FOR  Y=1  TO  Z 
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00630  D(Z/Y)=C(XCZ)-CY-1)"H)}!T(Z>Y-1)-(X(Z)-Y5SI-0!!T(Z,Y))/(2}{R(Z)£2) 

00640  NEXT  Y 

00650  NEXT  Z 

00660  FOR  Z=1  TO  J 

00670  FOR  Y=1  TO  Z 

00680  DCZ, Y)=-DCZ, Y)+BS2*ACZ, Y)+CX(Z)-CY-1)*H)*BCZ, Y)/R(Z) 

00690  NEXT  Y 
00700  NEXT  Z 
00710  FOR  Z=1T0J 
00720  FOR  Y=1T0Z 

00730  DCZ,Y)=CD(Z,Y)+BS2!{A(Z,Y)+RP(Z)“B(Z,Y)) 

00740  NEXT  Y 
00750  NEXT  Z 

00760  G(0)=0:GP(l)rl/Da,l):U=0 
00770  FOR  Z=2TOJ 
00780  W=0:V=0 
00790  U=U+GP(Z-1) 

00800  GCZ-1)=H«U 
00810  FOR  Y=1T0Z 
00820  W=W+0(Z,Y-1)S«GPCY-1) 

00830  V=V+BCZ#Y)S«GCY-1) 

00840  NEXT  Y 

00850  GPCZ)=C1-V/RCZ)-W)/DCZ,Z) 

00860  NEXT  .Z 
00870  FOR  Z=1T0J 
00880  U=0:V=0 


00890  FOR  Y=1T0Z 
00900  U=U+B(Z,Y)“GP(Y) 

00910  V=V+A(Z,Y)!*GP(Y) 

00920  NEXT  Y 

00930  JX(Z)=U 

00940  JR(Z)=VsiBS2 

00950  JO(Z)=l-dR(Z)-RP(Z)!\JX(Z) 

00970  NEXT  Z 
00990  FOR  Z=1T0J 

01010  GNCZ)=(C1+RPCZ)S2)!{C1-IXCZ))4WS2»IXCZ))s«JXCZ) 

01020  GACZ)=-15{MS2}«IXCZ)S2+C1+«PCZ)S2)“C1-IX(Z))S2 

01030  GBCZ)=CMS2-1-RP(Z)S2)5{JXCZ)S2 

01040  GC(Z)=Cl+JOCZ))S2 

01050  NEXT  Z 

01055  CLS 

01060  L PRINT* 'ENTER  THE  ANGLE  OF  ATTACK  A  CIN  RADIANS)”; : INPUT  A:LPRINT 
A 

01070  LPRINT' ’ENTER  THE  RADIAL  ANGLE  THETA  CIN  RADIANS)"; : INPUT  T:LPRIN 
T  T 

01071  LPRINT* 'ENTER  TEMP  CKELVIN)"; : INPUT  TE: LPRINT  TE 
01072  LPRINT"ENTER  PRESSURE  CN/M82)"; : INPUT  PE: LPRINT  PE 
01073  LPRINT" ENTER  DENSITY  (KG/MS3)"; : INPUT  PO: LPRINT  PO 
01080  FOR  Z=1T0J 

01090  CP(Z)=GOCZ)-2*'«GNCZ)5«SINCA)!{COS(A)«COSCT)+CGACZ)4GBCZ)«CCOS(T))S2 
-GCCZ)}«(SINCT))S2)5SCSINCA))S2 

01095  PR(Z)=l+200 . 9415*POs«TEJ{MS2}SCPCZ)/PE 

01100  NEXT  Z 
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01120  LPRINTTA8C1),,J,,;TABC10),,XCJ),,;TABC20),,CPCX),,;TABC32),,PRESSURE  RA 
TIOH 

01130  FOR  Z=1T0J 

01140  LPRINTTABa)Z;TABCl0)XCZ);TABC20)CPCZ);TABC32)PRCZ) 

01160  NEXT  Z 

01170  LPRINT"IF  YOU  DESIRE  TO  CHANGE  A  AND  THETA  ENTER  Y  ELSE  N";:INPU 
T  N$:LPRINT  N$ 

01180  IF  N$="Nn  THEN  END  ELSE  GOTO 105 5 
01200  XCZ)=XCZ>XCZ)«.001 
01210  RESUME 


Sample  data  Run  for  Mach  2.0  at  Angle  of  Attack  a  =  0  and  10  degrees 


SUPERSONIC  FLOW  PRESSURE  COEFFICIENTS  ABOUT  A  YAWED, 
SLENDER,  POINTED,  BODY  OF  REVOLUTION.  PROGRAM  BY 
J.R.SCH0N8ERGER  AM)  K.D.TILLOTSON.  BASED  ON 
CALCULATION  PROCEDURES  BY  PROF.  T.H.GAWAIN. 

NAVAL  POSTGRADUATE  SCHOOL,  MONTEREY  CA.  FALL  1980 

ENTER  XI  100 

ENTER  X2  260 

ENTER  X3  400 

ENTER  B1  .17453 

ENTER  B2  .0523599 

ENTER  COEFF  A-9.26731E-03 

ENTER  COEFF  B  8.36762 

ENTER  COEFF  C-402.881 

ENTER  RX2  33.856 

ENTER  FREE  STREAM  MACH  NBR  2 

ENTER  NBR  OF  SUBINTERVALS  20 

ENTER  LENGTH  OF  SUBINTERVAL  15 


J 

XCJ) 

R(J) 

RP(J) 

1 

21.5952 

3.80776 

.176324 

2 

43.1905 

7.61551 

.176324 

3 

64.7857 

11.4233 

.176324 

4 

86.3809 

15.231 

.176324 

5 

109.828 

20.1081 

.157449 

6 

129.756 

22.9529 

.129889 

7 

148.696 

25.2279 

.111218 

8 

166.989 

27.1293 

.0971743 

9 

184.81 

28.7579 

.0859282 

10 

202.262 

30.1736 

.0765366 

11 

219.414 

31.4156 

.0684513 

12 

236.311 

32.511 

.0613283 

13 

252.989 

33.4798 

.0549372 

14 

269.503 

34.354 

.0524078 

15 

286 

35.2186 

.0524078 

16 

302.498 

36.0832 

.0524078 

17 

318.996 

36.9478 

.0524078 

18 

335.493 

37.8124 

.0524078 

19 

351.991 

38.677 

.0524078 

ENTER  THE 

ANGLE  OF  ATTACK  A  (IN  RADIANS)  0 

ENTER  THE  RADIAL  ANGLE  THETA  (IN  RADIANS)  0 
ENTER  TEMP  (KELVIN)  288.15 
ENTER  PRESSURE  (N/MS2)  101325 
ENTER  DENSITY  (KG/MS3)  1.225 


J  X(J)  CP(X)  PRESSURE  RATIO 


1 

21.5952 

.114208 

1.31979 

2 

43.1905 

.114208 

1.31979 

3 

64.7857 

.114208 

1.31979 

4 

86.3809 

.114208 

1.31979 

5 

109.938 

.0989367 

1.27703 

6 

129.756 

.0651711 

1.18248 

7 

148.696 

.0455148 

1.12744 

8 

167.156 

.0324173 

1.09077 

9 

184.81 

.0229932 

1.06438 

10 

202.464 

.015794 

1.04422 

11 

219.414 

.0100775 

1.02822 

12 

236.311 

5.30874E-03 

1.01486 

13 

252.989 

1.29649E-03 

1.00363 

14 

269.503 

1.17926E-03 

1.0033 

15 

286 

3.38801E-03 

1.00949 

16 

302.498 

5.17499E-03 

1.01449 

17 

318.996 

6.62418E-03 

1.01855 

18 

335.493 

7.80804E-03 

1.02186 

19 

352.343 

8.76572E-03 

1.02454 

IF  YOU  DESIRE  TO  CHANGE  A  AND  THETA  ENTER  Y  ELSE  NY 
ENTER  THE  ANGLE  OF  ATTACK  A  (IN  RADIANS)  .0349 
ENTER  THE  RADIAL  ANGLE  THETA  (IN  RADIANS)  3.1416 
ENTER  TEMP  (KELVIN)  288.15 
ENTER  PRESSURE  (N/MS2)  101325 
ENTER  DENSITY  (KG/MS3)  1.225 
J  X(J)  CP(X)  PRESSURE  RATIO 


1 

21.5952 

.132234 

1.37026 

2 

43.1905 

.132234 

1.37026 

3 

64.7857 

.132234 

1.37026 

4 

86.3809 

.132234 

1.37026 

5 

109.938 

.117314 

1.32849 

6 

129.756 

.0821772 

1.2301 

7 

148.696 

.0616945 

1.17275 

8 

167.156 

.0482622 

1.13514 

9 

184.81 

.0377925 

1.10582 

10 

202.464 

.0303393 

1.08495 

11 

219.414 

.0236319 

1.06617 

12 

236.311 

.018247 

1.05109 

13 

252.989 

.0136563 

1.03824 

14 

269.503 

.0130074 

1.03642 

15 

286 

.0147623 

1.04134 

16 

302.498 

.0161278 

1.04516 

17 

318.996 

.0171903 

1.04813 

18 

335.493 

.0180226 

1.05046 

19 

352.343 

.0190157 

1.05325 

IF  YOU  DESIRE  TO  CHANGE  A  AND  THETA  ENTER  Y  ELSE  NN 


APPENDIX  B:  COMPUTER  PROGRAM  AND  SAMPLE  DATA  RUN 
FOR  CALCULATING  PROBE  PLENUM  PRESSURE 


Program 
00010  CLS 

00020  INPUT  "ENTER  PI"; PI 

00030  INPUT  "ENTER  P2";P2 

00040  INPUT  "ENTER  P3";P3 

00050  PRINT  "T=";T,"PP=";PP,"CASE=";CA 

00060  INPUT  "ENTER  YOUR  GUESS  FOR  PP";PP 

00070  A=P1-PP 

00080  B=P2-PP 

00090  C=P3-PP 

00100  IF  B<=0  AND  C<=0  THEN  200 

00110  IF  C<=0  THEN  300 

00120  IF  B<=0  THEN  400 

00130  PRINT'SOMETHING  WRONG"; : STOP 

00200  T=SQRC2,tABSCA))-25<SQRC25tABSCB))-SQRC2,sABSCC)) 

00210  CA=3 

00220  GOTO  50 

00300  T=SQRC2«ABSCA))+2«SQRC2stABSCB))-SQRC2»«ABSCC)) 
00310  CA=2 
00320  GOTO  50 

00400  T=SQRC25«ABSCA))-2«SQRC2s!ABSCB))+SQRC2}<ABSCC)) 

00410  CA»1 
00420  GOTO  50 


(Case  3) 
(Case  2) 
(Case  1) 


Sample  Data  Run  (Reproduced  from  CRT  display) 


ENTER  PI  ?1. 019237 
ENTER  P2  ?. 96525 
ENTER  P3  7.964412 
T=  0  PP=  0  CASE=  0 
ENTER  YOUR  GUESS  FOR  PP  7.977 
T=  .174618  PP=  .977  CASE=  3 
ENTER  YOUR  GUESS  FOR  PP  7.970 
T=  .0131546  PP=  .970  CASE=  3 
ENTER  YOUR  GUESS  FOR  PP  7.971 
T=-. 0186603  PP=  .971  CASE=  3 
ENTER  YOUR  GUESS  FOR  PP  7.9705 
1  T=-3 .07509E-03  PP=  .9705  CASE=  3 
ENTER  YOUR  GUESS  FOR  PP  7.9704 
T=  1.14933E-04  PP=  .9704  CASE=  3 
ENTER  YOUR  GUESS  FOR  PP  7.970404 
T=-1.11237E-05  PP=  .970404  CASE=  3 


APPENDIX  C:  SEQUENTIAL  LISTING  OF  TESTING 


As  stated  in  Chapter  3,  Appendix  C  is  a  sequential  listing  of  testing 
and  corrective  action  taken  between  testing  attempts,  and  a  discussion  of 
the  three  most  likely  reasons  for  the  failure  of  the  wind  tunnel  to  start 
at  Mach  1.4.  Testing  attempts,  and  the  corrective  actions  taken  between 
testing  attempts,  were  as  follows: 

(1)  Test  Number  One  was  conducted  with  the  Mach  1.4  nozzle  blocks 
installed  and  the  model  mounted  in  the  wind  tunnel  test  section.  Super¬ 
sonic  flow  could  not  be  established  in  the  test  section.  Test  section 
blockage  was  considered  to  be  the  most  probable  cause.  A  higher  percent¬ 
age,  28%,  of  test  section  blockage  is  allowed  at  Mach  2;  therefore,  the 
decision  was  made  to  complete  Mach  2  testing  prior  to  reducing  model  and 
mounting  system  frontal  area.  The  Mach  1.4  nozzle  blocks  were'  removed. 

The  wind  tunnel  plenum  chamber  screen  was  discovered  to  be  torn,  and  was 
hanging  loose  in  the  plenum  chamber.  The  two  functions  of  the  screen  are: 
(a)  to  obtain  a  flat  velocity  profile  (flat  stagnation  pressure  profile) , 
and  (b)  to  generate  small-scale  turbulence  which  dissipates  prior  to  ar¬ 
rival  at  the  test  section.  Sufficient  manpower  to  replace  the  plenum  cham¬ 
ber  screen  was  not  available.  Therefore,  the  screen  was  removed. 

(2)  The  Mach  2.0  nozzle  blocks  were  installed.  To  insure  that  the 
flow  in  the  wind  tunnel  test  section  was  uniform,  an  additional  pressure 
tap  was  installed  in  the  wind  tunnel  wall  opposite  one  of  the  existing 
pressure  taps.  Testing- showed  no  significant  pressure  difference  (approxi¬ 
mately  0.5  percent)  between  the  two  pressure  taps.  Therefore,  the  test 
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section  velocity  profile  was  inferred  to  be  reasonably  uniform.  Twenty- 
three  Mach  2.0  wind  tunnel  tests  were  conducted  with  supersonic  flow 
achieved  in  the  test  section  (25  to  30  PSIG  plenum  chamber  pressure). 

(3)  The  Mach  2.0  nozzle  blocks  were  removed,  and  the  Mach  1.4  noz¬ 
zle  blocks  were  reinstalled.  The  model  was  mounted  in  the  wind  tunnel, 
and  testing  was  attempted.  Supersonic  flow  in  the  test  section  was  not 
achieved. 

(4)  The  model  frontal  area  was  reduced  by  14  percent  (yielding  5.4 
percent  test  section  blockage)  by  machining  0.050  inch  from  the  top  and 
bottom  of  the  missile  nose  ogive.  Supersonic  flow  still  could  not  be 
achieved. 

(5)  Model  frontal  area  then  was  reduced  by  an  additional  30  percent 
(yielding  3.5  percent  test  section  blockage)  by  machining  an  additional 
0.060  inch  from  the  top  and  bottom  of  the  missile  nose  ogive  and  0.015 
inch  from  the  top  and  bottom  of  the  double  edge  airfoil.  After  the  reduc¬ 
tion  in  blockage,  supersonic  flow  could  not  be  achieved. 

(6)  Four  boundary- layer  bleed  valves  were  installed  in  both  wind 
tunnel  walls,  just  forward  of  the  nozzle  second  throat,  to  decrease  the 
required  mass  flow  rate  through  the  second  throat.  A  complete  description 
of  the  bleed  valve  specifications  and  function  can  be  found  in  Babel  [8] . 
Once  again,  supersonic  flow  could  not  be  achieved. 

(7)  The  bleed  valves  and  the  model  were  removed  from  the  wind  tun¬ 
nel.  Attempts  to  obtain  supersonic  flow  in  the  wind  tunnel  test  section 
failed  with  no  wind  tunnel  blockage. 

(8)  A  fine-mesh  screen  was  installed  in  the  wind  tunnel  just  for¬ 
ward  of  the  nozzle  first  throat.  With  the  screen,  but  without  the  model. 
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supersonic  flow  could  not  be  achieved.  Figure  A-l  is  a  photograph  of  the 
shock  that  was  observed  in  the  wind  tunnel  test  section  between  25  and  65 
PSIG  plenum  chamber  pressure. 

The  three  most  probable  causes  of  the  failure  to  obtain  supersonic 
flow  were:  (a)  that  the  model  was  too  large,  resulting  in  choked  flow;  (b) 
that  the  wind  tunnel  plenum  chamber  screen  was  missing,  causing  an  uneven 
velocity  profile;  and  (c)  that  the  Mach  1.4  nozzle  block  second  throat  was 
too  small,  resulting  in  choked  flow.  These  factors,  of  course,  can  operate 
in  combination  as  well  as  singly. 

Model  size  was  ruled  out  as  a  cause  of  failure  due  to  two  considera¬ 
tions.  First,  after  reduction  of  the  model's  frontal  area,  the  actual  wind 
tunnel  blockage,  3.5  percent,  was  well  below  the  6.4  percent  theoretical 
allowed  wind  tunnel  blockage.  Second,  supersonic  flow  could  not  be  achieved 
even  without  the  model  in  the  wind  tunnel.  The  wind  tunnel  plenum  chamber 
screen  and  the  nozzle  block  second  throat,  however,  remained  as  probable 
causes  of  testing  failure. 

The  inside  dimensions  of  the  Mach  1.4  nozzle  blocks  were  measured, 
and  the  areas  of  the  first  and  second  throats  were  14.36  and  15.61  square 
inches,  respectively.  The  ratio  of  second  throat  area  to  first  throat  area 
was  thus  1.09.  Using  the  calculation  method  presented  in  Liepman  and  Roshko 
[6],  the  theoretical  minimum  area  ratio  allowable  is  1.04.  Therefore,  the 
Mach  1.4  nozzle  block  design  allows  a  safety  factor  of  approximately  five 
percent  for  the  boundary  layer  in  the  second  throat. 

Since  there  is  less  constriction  in  the  Mach  1.4  nozzle  than  in  the 
Mach  2.0  nozzle,  an  uneven  velocity  profile  would  affect  more  adversely 
the  flow  in  the  Mach  1.4  nozzle.  The  adverse  effect  of  the  uneven  velocity 
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00470  w=w+acz,y;3'fpcy; 

00480  NEXT  Y 
00490  IXCZ)=W 
00500  W=0 
00510  NEXT  Z 
00520  FOR  Z=1T0J 
00540  NEXT 
00550  FOR  Z=1T0J 

00560  GOCZ)=Cl-a+RPCZ)S2)5«(l-IXCZ))S2+<M»«IXCZ))62) 
00565  NEXT  Z 
00610  FOR  Z=1  TO  J 
00620  FOR  Y=1  TO  Z 

65 


END 
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